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Abstract 

> 

' A three-dimensional reduction of the two-particle Bethe-Salpeter equation is 

proposed. The proposed reduction is in the framework of light-front dynamics. 

G\ ' 

It yields auxiliary quantities for the transition matrix and the bound state. 
The arising effective interaction can be perturbatively expanded according to 
the number of particles exchanged at a given light-front time. An example 
suggests that the convergence of the expansion is rapid. This result is par- 
^ ' ticular for light-front dynamics. The covariant results of the Bethe-Salpeter 



equation can be recovered from the corresponding auxiliary three-dimensional 
ones. The technical procedure is developed for a two-boson case; the idea for 
an extension to fermions is given. The technical procedure appears quite 
practicable, possibly allowing one to go beyond the ladder approximation 
for the solution of the Bethe-Salpeter equation. The relation between the 
three-dimensional light-front reduction of the field-theoretic Bethe-Salpeter 
equation and a corresponding quantum-mechanical description is discussed. 
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I. INTRODUCTION 



In relativistic field theory the Bethe-Salpeter equation (BSE) |T| describes two-particle 
systems in interaction. The inhomogeneous BSE 

T — V + VG T (1) 

yields the transiton matrix T of two-particle scattering. In Eq.(|l]) Gq is the free resolvent 
which propagates two non-interacting particles, i.e., 

Go = - - , (2) 

k\ — m\ + iok\ — 777J + io 

fcf denoting the momentum operator of particle % with mass m,, the hat on the variable 
emphasizing its operator character. The driving term V stands for the complete interaction, 
irreducible with respect to two-particle propagation; it also includes self-energy corrections, 
i.e., it may contain disconnected pieces. If the dynamics allows for a two-particle bound 
state with total four- momentum Kb , K B = M B , the vertex |T) at the bound-state pole 
is solution of the homogeneous BSE 

|r> = vg |r> (3) 

with the relation 

1*) = Go m (4) 

to the bound state |^). Eqs. (|1|) and (|3]) do not determine the bound state |\&) in full; 
the normalization condition has to be added. The two-particle total four-momentum K is 
conserved, i.e., all operators O a of Eqs. ([I]) and (^), O a = T,G or V, and the states |^/) 
and |T) carry a four-dimensional 5— function in momentum space, i.e., 

(K'\O a \K) = S(K'-K)O a (K), (5) 
(K'\ *) = 5(K , -K B )\V B ), (6) 



(K'\ T) = 8{K' - K B ) \Y b ) , (7) 

the reduced quantities depending parametrically on K, even if not spelled out explicitly 
for the states |r#) and \^b)- The reduced states \^b) and \Tb) belong to and the oper- 
ators O a (K) act in a Hilbert space characterized by a four-dimensional momentum or 
coordinate x M . They satisfy the Eqs. ([T]) and ([J) in a corresponding fashion. 

The inhomogeneous and homogeneous BSEs ([]]) and (f|) are general and exact formula- 
tions for the scattering amplitude and bound state. However, for any realistic field theory 
solution of the BSE constitutes a difficult calculational task which has not been tackled in 
full. In practical calculations, the driving term V(K) has to be truncated to low orders 
of particle exchange. In Euclidean space, the fermion case has only been solved in ladder 
approximation 0, i.e., with single particle exchange for the driving term while the boson 
case has only been solved in ladder and crossed ladder approximation ||. However, the 
step from the Euclidean-space to Minkowski-space solutions requires a complicated analytic 
continuation Q. Direct solutions in Minkowski space are just now becoming available ||. 

In the light of the great calculational difficulties, three-dimensional reductions of the BSE 
are still of high physics interest. The conceptual sacrifices generated by the reduction can 
possibly be outweighed by the gain in technical ease: One hopes to be able to include physical 
phenomena which the four- dimensional BSE with a highly truncated interaction is unable to 
account for. For example, the three-dimensional Gross approach allows only one particle 
to propagate off- mass-shell, but it appears to go beyond the ladder approximation of BSE by 
single particle exchange and to include crossed exchanges implicitly; it manfestly preserves 
covariance. Other reduction schemes give up covariance, which then must be recovered 
through complicated correction schemes. An equal-time projection scheme has also been 
explored for the pion-nucleon system which fulfills requirements of covariance and discrete 
Poincare symmetries 0. The papers by Fuda [Q report on the comparision of one-meson 
exchange models in ladder approximation on both light-front and instant-form dynamics, 
without emphasis to the underlying field-theoretic framework. 
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The purpose of this paper is two-fold: 

i) First, the paper attempts to find a three-dimensional equation for auxiliary quantities 
from which the full covariant solution of the BSE in the ladder or any other approximation 
can be obtained with ease. This is a technical objective with solutions well-known in the 
framework of instant-form dynamics. Here the advantages of light-front dynamics are to be 
explored. 

ii) Second, the paper tries to illuminate the connection to a quantum-mechanical de- 
scription of the two-particle system whose dynamic input is related to the underlying field 
theory. 

Sect. II motivates our novel choice for three-dimensional auxiliary quantities from which 
the covariant solutions of the BSE are obtained. It motivates light-front dynamics as our 
choice for a dynamical framework. Sect. Ill gives our theoretical apparatus in full. Sect. IV 
tests the potential of the method in the example of a two-boson bound state. We perform 
numerical calculations for the two-boson bound state including up to four-particle intermedi- 
ate states in lowest order and compare to the solutions of the four-dimensional BSE equation 
in the ladder approximation. Sect.V sketches the generalization of our theoretical apparatus 
to fermions. Sect. VI discusses the connection with light-front quantum mechanics. Our 
conclusions are summarized in Sect. VII. 

II. CHOICE OF TWO-PARTICLE AUXILIARY FREE RESOLVENT G Q {K). 

It is well known, from the work of Ref. f9j, that the transition matrix T(k) and the 
bound state \^b) of the covariant BSE can be obtained with the help of a convenient 
auxiliary resolvent Gq(K), still to be chosen. That is, we have 



T(K) = W{K) + W{K)G {K)T(K), 



(8) 




(9) 



* B > = G (K) \T B ) , 



(10) 
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provided the driving term V(K) is changed to W(K) according to 
W(K) = V(K) + V(K)[G (K) - G (K)]W(K) 



(11) 



Eqs. (0) and (O) do not determine the bound state \^ B ) m hill; the normalization condition 

G (K)-i - G (K B yi V(K) - V{K B ) \ = i 



lim ( ^ B 

R2^Kl \ 



K 2 -K 2 B 



K 2 -K 2 B 



has to be added. It involves the original driving term V(K) JTIJ. The choice of Gq(K) is 
hoped to be sufficently clever that the integral equation ( P]) does not have to be solved in 
full, but that a few terms of the infinite series 



W(K) = V(K) ]T [(G (K) - G {K)) V(K)] n 

n=0 

W(K) = V(K) + V(K) (G (K) - G (K)) V{K) + .. 



(13) 



suffice. The auxiliary resolvent Gq(K) remains a four- dimensional one, but its choice may 
sacrifice the covariance which the resolvent G (K) possesses. 

The dynamics of the interacting two-particle system can be fully described by its prop- 
agation between hyperplanes, the hyperplanes x°= const, in instant-form dynamics, the 



hyperplanes x 



x + x 



const, in light-front dynamics |TT|. In contrast, the free resol 



vent of the BSE depends on the individual times x® or on the individual light-front times 
xf. 

The free resolvent in instant-form coordinates ki = (h®, k^) 



~'0 '0 

Xi x 2 



G(\ 



~o o 
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(27T) 

e 



^ / dk^dK x 

7T ) J 



lri *a "l T *2/ 



(k^f-k -m 2 + io\ ( {K° — k^) 2 — (K — ki) 2 — m 2 , + io 



(14) 



in fact only its dependence on individual times x® is made explicit reduces for 

propagation between the hyperplanes x° and x'° to 

dK° 



x'°x'° 



Go 



x°x° 



iK u x'"-x u 



27T 



dk^dkUk? G (AT) k" 



(15) 
(16) 



In Eg. (|15|) the notation 
k? G (K) = _ 1 ^(fcf-fc?) 



2tt 



(A;?) 2 - fci - ml + io) ((K° - k\) 2 - (K - h) 2 - m 2 + io^j 



is introduced, as well as the abbreviation 



,G (K)\ : = J dk'?dkl (k?\ G (K) (18) 

■ (19) 



2k° lon 2k° 2on \ (Jfo - A;? on - fc 2 ° on + io) [K° + k° lon + k\ on - *o | 

The matrix element {kf\ Gq(K) of Eg. flPTD , in which only the dependence on the "dy- 
namic" variable k® is made explicit remains an operator with respect to the "kinematic" 
variables hi, the operator character being carried by the operators k\ on = y k i + m 2 . The 
basis states for these kinematic variables are defined by (Si ki) = exp (iki ■ Xi) and are 
eigenfunctions of the momentum operator k and the free energy operator k° on . The states 
k^i form an orthogonal and complete basis. 

In Eg. fll8|) , the vertical bar |o indicates that the dependence on is integrated out. 
The bar on the left of the resolvent represents integration on k® in the bra-state, the bar 
on the right in the ket-state; we shall encounter resolvents in which integration on is 
performed only on one side, the bar |o being placed on that side alone. The resulting operator 
|oGo(-K")|o is three-dimensional and depends only on the kinematic variables k\. It is a global 
propagator, since it mediates between hyperplanes according to Eg. (pT6|) , not allowing for 
individual time differences between the two particles, it is not explicitly covariant. In instant- 
form dynamics, the global propagator | Cro(^0|o still allows for particle and antiparticle 
propagation. This is considered to be a technical disadvantage. 

The free resolvent in light- front coordinates ki = (kj~ := k° — k 3 , kf := k° + k 3 , k±) 



-'+,_/+ 

1 2 



Go, 
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only its dependence on the individual light-front "times" xf is made explicit 

for propagation between the hyperplanes x + and x' + , to 

dK~ 



xx 



\K-(x'+- 



2tt 
dK~ 



dkffdki (k[- G {K) fcT/) , 



e-i K ~^-^)\G (K)\ ■ 



reduces, 

(21) 
(22) 



In Eq. 



the notation 



k'r G (K) hi 



5 (k'f - kl 



2tt 



k2±+m2—io 
K+-kt 



(23) 



is introduced with the abbreviation, 



\G (K)\: = J dk[~dk^ (k[-\G (K)\kl 

i6(K+ - %t)9{ki 
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ki(K + — ki ) yK — k lon — k 2on + 
:= 9o(K) , 



(24) 
(25) 
(26) 



where K + > can be chosen without any loss of generality. The matrix element 
(k'{~ Gq(K) ki \ of Eq.(p3|), in which only the dependence on the "dynamic" variable fcf is 
made explicit, still remains an operator with respect to the "kinematic" variables (kf, ki±), 



and k 



2on 



(K±-k 1± )' 2 +m 2 



variables are defined by 



\X^ Xi 



K+-kQ 



k i ki 



-. The basis states for the kinematical light-front 



(27) 



and are eigenfunctions of the momentum operators {kf ,ki±) and the free energy operator 



kxr.y, . The states 



k + k±) form an orthonormal and complete basis, e.g., 



dk + d 2 k 



(x' x' ± k + k±j (k + k± x x±\ = 5(x' — x )S(x' ± — x±) 



(28) 



2(2vr) 3 

In Eq. (El) the vertical bar | indicates that the dependence on fcj" is integrated out. The 
bar on the left of the resolvent represents integration on fcf in the bra-state, the bar on the 



right in the ket-state. We shall encounter resolvents in which integration on k± is done on one 
side alone, the bar | being placed only on that side. The operator Qq{K) is three-dimensional 
and it depends on the kinematic variables (kf , k\±) only. It is a global propagator, since it 
mediates between hyperplanes according to Eq. (|2!^) , not allowing for individual light-front 
time differences between the two particles. It does not possesses explicit covariance but is 
still covariant under light-front boosts. In light-front dynamics, the global propagator go(K) 
only allows particle propagation, no antiparticle propagation, due to the choice of K + > 0. 
This is the advantage of light-front dynamics, with which we work from now on. 

The auxiliary four- dimensional resolvent Gq(K), introduced in Eqs. (|8|)-(|T3l) has to be 
chosen next. We require for G (K): 

G (K)\ = G (K)\ , (29) 
\G (K) = \G (K) , (30) 
\G (K)\ = \G (K)\ , (31) 

and define a three-dimensional transition matrix t(K) through 

| [G (K) + G (K)T(K)G (K)] | = g (K) + g (K)t(K)g (K) . (32) 

In Eqs. (^9"|) -(|32l) the abbreviation | for integrating out the fcf dependence of operators is 
used. The conditions fl29D -(|32"D are a rather mixed bag. The conditions (31) and ( |32"1) are 



physical ones: They require that the global-propagator form of Gq(K) be the same as for the 
exact free resolvent Go(K) and that the full resolvent of BSE Gq{K) + Gq{K)T{K)Gq{K) 
can be obtained from |Go(-^')| and the three-dimensional t(K). However, the two conditions 



(pl|) and (i32|) do not determine Gq(K) in full. Our choice is 

G (K) := G (K)\g^(K)\G (K) , (33) 

though G (K) = 5 (k'{~ — 9o{K)S (ki — (and obvious variants of it) seems to be 
a legitimate alternative. However, if we demand that the kernel of the integral equation 
for the auxiliary transition matrix, t(K), represents light-front propagation in higher Fock- 



states, then the choice is unique. The conditions (|29D and (|30D introduce the additional 



convenience that the auxiliary resolvent be as close as possible to the exact free one. The 
auxiliary quantities are computed in Appendix A. 



III. CALCULATIONAL PROCEDURE 

Our calculational procedure amounts to solving three-dimensional integral equations, 
whose solutions then yield the covariant results of the BSE by quadrature. 

The four-dimensional transition matrix T(K) is obtained from the three-dimensional 
auxiliary one t(K), defined by Eq. fl32|) , through 



t(K) = g,{K)- l \G,{K)T(K)G,(K)\g,(K)- 1 , (34) 

by first iterating the integral equation (§) once, 

T(K) = W{K) + W(K) [G {K) + G (K)T(K)G (K)] W{K) , 

and then making use of our choice, Eq . (|33|) , for Gq(K) and the result Eq.(^). The relation 
between the T(K) and the auxiliary t(K) is 

T{K) = W{K) + W(K)G (K)\ [g^K)' 1 + t{K)] \G (K)W(K) . (35) 

The auxiliary transition matrix t(K) itself is obtained by the three-dimensional integral 
equation 

t(K) = w(K) + w(K)g (K)t(K) , (36) 

in which the driving term w(K) is derived from the modified four-dimensional interaction 
W(K) of Eq.(p]) according to 

w(K) := goiKy'lGoWW^GoWlgoiKy 1 . (37) 

There is an integral equation for w(K) as there is for W(K), but we do not give it here. 
We hope that, through our choice (^) for Gq(K), a few terms of the expansion Eq. (|TT|) , 
of W(K) in powers of V(K) will dynamically suffice to yield the full result of BSE with 
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satisfactory accuracy. The numerical example of Sect. IV where rapid convergence of w(K) 
is seen, demonstrates the validity of this expectation. 

If the transition matrix T{K) of the BSE has a bound-state pole at total four momentum 
K B , K b = M|, the auxiliary three-dimensional transition matrix t(K) also has a bound- 
state pole at exactly the same K B , according to Eq. fl3~4]) , with the residue \jg) being the 
solution of the homogeneous three-dimensional equation 

\ lB ) = w{K B )g {K B ) \ lB ) , (38) 



corresponding to the inhomogeneous one, Eq . (p6|) . From \j B ), the residue \T B ) of BSE can 
be recovered according to Eq. (|35|) 



\T B ) = W(K B )G (K B )\\ lB ) (39) 

as well as the bound state \^? B ) of BSE, i.e. 

\* B ) = G (K B )W(K B )G (K B )\ \ lB ) , (40) 

\V B ) = [l + (G (K B ) - G (^b)|^o(^)" 1 |G'o(^)) W(K b )] G (K b )\ \j b ) . (41) 



For the form Eq.(4"ip of the bound state, the condition Eq. (53) \ j B ) — w(K B )g (K B ) \ j B ) = 
is used. The step from the three-dimensional residue to the four- dimensional bound 
state \^ B ) appears predominantly a kinematic one, effected by the operator Gq(K b )\ . Only 
the second term in Eq . (|4lD depends on the interaction, and it is expected to be a small 
correction. 

The four-dimensional bound state \^f B ) is related to the auxiliary three-dimensional \4> B ), 
defined by 

\<Pb) := g (K B ) \ lB ) (42) 

and satisfying 

\<Pb) = g (K B )w(K B ) \<j> B ) , (43) 
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in an obvious way by 



J dki(k^ B ) = \4> B ) . 



(44) 



The result Eq. (j44|) follows immediately from Eq. (j4lr) . The auxiliary bound-state wave- 
function \<p B ) is the projection of the bound-state \^b) of BSE to equal light-front individual 
times xf = x + , taken on the hyperplane x + = 0. 

The bound-state \^b) of BSE and its three-dimensional auxiliary version \4> B ) still have 
to be normalized. If the dependence on K of the original interaction V(K) is weak, i.e., 
(V(K) — V(K B ))/(K 2 — K B ) — and if furthermore the interaction-dependent term in the 
step from \<p B ) to \^b) according to Eq.(^TJ) is small, i.e., \^ B ) — Go(KB)\go(KB)~ 1 \4>b), 
then 



lim ( 


'^B 


Go(K)- 1 


-G (Kb)- 


1 


V(K) 


- V(K B ) 


K 2 


- K 2 

t^B 




K 2 


- K 2 


lim ( 


^B 


Go(K)- 1 


-G (Kb)- 


1 


^B^j ^ 




K 2 


- K 2 

t^B 






lim ( 


^4>B 


go(K)-' - 


go(KB)- 1 




b\ = 1 • 




K 2 - 


- K 2 


4 





B 



(45) 

For any further applications, i.e., for predicting physical observables, we now have two 
equally valid options. We may either work with covariant operators using the bound state 
\^b) and/or the transition matrix T(K) of the BSE or we may derive effective operators 
suited for the context of the auxiliary three-dimensional bound state \4>b) and/or the aux- 
iliary three-dimensional transition matrix t(K). We give an example of each of the possible 
strategies: 

We use the eletroweak current J^{Q) as example and assume that it connects an initial 
bound state \^m) to a final one \^ B f) in an elastic process. We take J^{Q) to be the current 
appropriate for the hadronic field theory with four- momentum transfer Q = Kb/ — K B i- 
The matrix element for describing the process (^Bf \ <J^{Q) \^Bi) can be obtained from the 
three-dimensional bound state \^ B ) by 



(* B/ | J»(K Bf - K Bi ) \*Bi) = (<j> Bf \f{K Bf , K Bi ) \<j> Bi ) 
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(46) 



with the effective current in three-dimensional space 

f{K f , Ki) := g Q {K f y l \G Q {K f ) [l + W{K f ) (G (K f ) - G Q {K f )\g Q {K f )- l \G Q (K } ) 
J»(K f - Ki) \l + (G (Ki) - GoiK^golK^lGoiK)) W(K)} G^K^g^K)- 1 . (47) 



For the relation between the bound states \^b) and \4>b), Eq.(f4"T|) is used, which separates 
the kinematic and dynamic, i.e., interaction-dependent, steps in that relation from each 
other. The bound state has to be calculated for the initial and final four-momenta K Bi 
and K B f- The effective current j^iK^ Ki) is predominantly derived kinematically from the 
covariant one through g (Kf)~ 1 \G (Kf)J" ti (Kf — K i )Go(K i )\g (K i )~ 1 but it also depends on 
the interaction W(K) of Eq. flTTD . If W(K) is not computed in full, but only expanded up 
to a certain order in the original interaction V(K) of the BSE, the effective current should 
be expanded consistently up that order. 

IV. A NUMERICAL TEST CASE 

We use the bound state of a schematic two-boson system as a test case of the power of 
the suggested numerical technique. The employed interaction Lagrangian is 

£/ = #s0i</>i0- + gs4>k<t>20, (48) 

where the bosons with fields 0i and 02 have masses mi and m.2, which we take to be equal, 
mi = m.2 = m, and the exchanged boson with field a has mass fi. The coupling constant is 

Using standard techniques in Euclidean space, the homogeneous BSE is solved for the 
bound-state vertex \T B ) in the ladder approximation, i.e., 

/u |r \ = • 2 f (fci \^b) ( , q) 

Ws J (2vr) 4 ((k{ - h) 2 -fi 2 + is) (kf - m 2 + ie) ((K B - h) 2 - m 2 + is) ' 1 ] 

The solution is calculated in the two-particle cm. system, i.e., for Kb = (Mb,0), and 
for the ratio of masses /i/m = 0.5. Requiring the bound state mass to have a particular 
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value Mb fixes the coupling constant g$- The four-dimensional bound-state vertex (ki \Tb) 
depends on all Euclidean four components of the momentum k\ of boson 1. The exact four- 
dimensional bound state is obtained according to Eq.(|H]). However, the representation of 
the vertex and bound state in terms of Minkowski momenta is difficult. We do not attempt 
it. 

In contrast, the four-dimensional bound-state obtained by the numerical technique sug- 
gested in Sect. Ill is available in Minkowski space. We calculate it only approximately by 
using for the driving term w(Kb) of the auxiliary three-dimensional equation Eq.(^), an ex- 
pansion in orders of the interaction V(K) of BSE in Eqs. ( |13D and (p7|), i.e., in powers of the 
coupling constant gs of the interaction Lagrangian (f48|). We use the approximation up to the 
second and fourth powers of gs, i.e., w(K B ) — w^ 2 '(Kb) and w(K B ) — w^ 2 \K B ) +w^(K b ). 
In a time-ordered view, the BSE allows for an exchange of an infinite number of a bosons 
in stretched configurations. In contrast, the approximative w^ 2 \Kb) allows only for one 
exchange (Fig.la), while w^(Kb) allows for two (Fig.lb). The analytic forms of w^(Kb) 
and w^(K B ) are given in Appendices B and C. The explicit forms of the homogeneous 
integral equation for 17b), Eq.(^), for the above approximations in the driving term are 
given in Appendix D. In order to make a comparision with the exact bound state we study 
the projected forms of the bound states, i.e., 

fe X act{\IWs) ■ = J dKdk+ih |tf fl > (50) 

= 2 J dk%dkl(k%k 1± kl \G (K B )\T B ) . 

/ig(V§J = / dkf(kfk 1± |<Mi;) (51) 
= Jdkf(kfk 1± \g (KB)\lB) { :J p . 

The superscripts (n) in Eq. fl5"T|) indicate the power of the coupling constant up to which the 
approximation is carried, i.e., w(Kb) — Ya=2 w (Kb)- The comparision between exact and 
approximate results is carried out on two levels: 
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In Fig. 2 the relation between gg and Mb is tested for // = 0.5m against the four- 
dimensional results. Whereas the exact relation is already satisfatory reproduced by the 
approximation based on w^ 2 \K B ), the approximation based on w^ 2 \Kb) + w^'(K B ) im- 
proves the agreement. 



In Figs. 3 and 4, the projected bound-states f(yk 2 ± ) are compared for two cases. In 
the first case Mb = 0, i.e., the binding is very strong. It is of the order of the masses of 
the interacting particles as encountered in quark systems. In the other case Mb = 1.98m, 
i.e., the binding is very weak. It is only 2% of the masses of the interacting particles, as 
encountered in nuclear systems. In both cases the approximation based on w^(Kb) is 
already quite accurate. The improvement due to the inclusion of w^(K B ) is particularly 
noticeable for the case of strong binding. 

The fact that a low-order approximation of w^ h \Kb) works surprisingly well is a virtue of 
light-front dynamics. It is well known that the analogous approximation scheme in instant- 
form dynamics has much poorer convergence properties with respect to the number of ex- 
changed a bosons [|T2]]. 

V. EXTENSION TO FERMIONS 

The free resolvent which propagates two fermions disconnectedly contains self-energy 
corrections as in the case of bosons. They are usually left out of the ladder approximation 
of interaction. The two-fermion free resolvent then takes the form which we immediately 
rewrite conveniently as 



G F = ?^5H^ (52) 



-iF _ $1 + m l + m 2 

k\ — m\ k\ — m 2 , 

Gl = AGq + (} lon + mx) (hon + m 2 ) G , (53) 



where k lon = fclJ ^+ m i anc l k 2on = k2J ^ m ' 2 . In Eq. (p3[) G is the covariant propagator the 

k 1 k 2 



paper has worked with in the conceptual development untill now. Furthermore, Eq.(|53"| 
is the definition of AG^ which contains - except for the particular spin-dependent opera 
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tors (jfc lon + mij and {fc 2 on + that commute with Go - all particular divergences and 
subtleties connected with the fermion motion. The operator 



A F = 7l hon + m 2 hon + m l it , 7l it 

" ~ 2fc+ fc2 _ m 2 + kl-ml 2k+ ^ 2kf 2k+ ' ' 



carries the instantaneous part of the fermion propagators in light-front time. Its is singular 
under integration. We therefore suggest the following strategy for fermions: We apply 
the reduction to an auxiliary resolvent G twice, using the apparatus of Sects. I and II. The 
operator dependence on the total two-fermion four momentum K is factored out as there. 
All operators become then parametrically dependent on K. 

In the first step, the two-fermion resolvent (ft lm + md (^ 2 on + m s) G (K) is introduced 
instead of Gq(K). We use formulae fl8|)-(|l3|) to do this. All the physics of anomalous two- 



fermion propagation is contained in the new effective interaction W(K) of Eq.flTTp. Thus, 
one arrives at a new BSE, corresponding to Eq.(|]) after reduction with respect to K, with 
the four dimensional resolvent (f> 1(m + mij f$ 2o n + m2 ) ^oC^O an d the new interaction. The 
resulting two-fermion equation is now solved with the technique as developed for two bosons. 
This is possible due to the fact that the spin-dependent operator (ft lon + mij (^ 2 on + 
also commutes with the auxiliary one G (K) = G (-^')|fl , o(-^')~ 1 |G ! o(-ft'), i.e., 

{hon + m,) (h on + ma) G (K) = G (K) (} lon + m x ) (fe on + m 2 ) . (55) 

This idea will not be further developed in this paper, but indicates that the scope of the 
method extends beyond the two-boson system. 

VI. RELATION TO LIGHT-FRONT QUANTUM MECHANICS 

Sects. I-IV used the notion of a bound state, but scattering states were not introduced. 
The later could have been introduced in the BSE as well as in the auxiliary three- 
dimensional equation ( |36| ) for t(K) with the global propagator go(K). Given an initial two- 
particle plane-wave state kfkij_K on ) with total momentum i^ on and light-front "energy" 
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K _ _ kl ± +m\ (j? x _fc 1± )2 +m 2 
on I K+-k+ ' 



one may define the corresponding three-dimensional scattering 



state 



4>( + \kiki±K on ) ) with outgoing light-front boundary conditions as the solution of a 



(56) 



standard Lippman-Schwinger type of equation, i.e., 



-\kfk 1± K 



kfk 1± K on ) + g Q (K on )w(K on ) ^ + \kfk 1± K, 



with four- momentum K on = (K on , K + , K±). The relation to the auxiliary transition opera- 
tor t(K) is obvious, 



t{K on ) \kfk 1L K on ) = w(K on ) \(P i+ \kfk 1± K on )) . 



(57) 



Furthermore, it satisfies the homogeneous equation 



(Kon)- 1 - w(K on )} \<P {+ \kfk l± K on )) = 



9o( 



in the same way as the auxiliary bound state \4>b) of Eq. (|38|) does, i.e., 



(58) 



g (K B )- 1 - w{K B )] |<M = 



(59) 



Eqs. (^) and (^) formally look like the eigenvalue equations of quantum mechanics with 
the only difference being that the two-particle interaction w(K) depends on the eigenvalue. 
Untill now the relationship to quantum mechanics has indeed been entirely formal. The 
states \4>b) and (kf ki^K on )^ and the corresponding transition matrix have significance 
only as quantities from which the solutions of the BSE can be obtained with comparative 
ease. On the other hand, at this stage a quantum-mechanical description of the two-particle 
system can be given which corresponds dynamically to the underlying field-theoretic one, 
though it is by no means equivalent to it. Quantum- mechanical two-particle states \<p) are 
required to satisfy the eigenvalue equation for the squared mass operator, 



M* + v(K + ,K ± ) \ip)=M^\ip) 



(60) 



where the squared free-mass operator is 



Ml 



-2 

k 1± + mj 



+ 



x 



- h±) 2 + m 2 2 
1 — x 



(61) 
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and x = j^. The states are elements of a Hilbert-space spanned by the free-particle on-mass- 
shell basis states. Boundary conditions must be imposed on the solutions of Eq. flBDD in order 
to make them acceptable. Bound-state and scattering state solutions to the mass squared 
operator equation exist and are orthonormalized. The ort honor malizat ion for scattering 
states is of the 5-function type. The states have a probability interpretation. The quantum 
mechanical bound-state normalization is 

(<Pb\<Pb) = 1 • (62) 

The two-particle potential v(K + , K±) is independent of the eigenvalue K~, the eigen- 
value Kg to be calculated for the bound state and the eigenvalue K~ n prescribed for the 
scattering states; the potential is hermitian it is instantaneus in light-front time; it conserves 
the kinematic components (K + , K±) of the total two-particle four-momentum K. In quan- 
tum mechanics v(K + , K±) may be parametrized by fitting it to observables. If contact is 
attempted to a corresponding field theory a standard form of identification is 



(k?k' 1± \v(K+,K ± )\ktk 1± ) := iJ y+^rz fc , +) (K + k'i± HK V )\ k + k x )^ ] ^^- fc + (63) 
with 

K v = (±K' m + ±K m ,K + ,K ± ) . (64) 

The defined relativistic quantum mechanical interaction is cooked in the framework of light- 
front dynamics. The value K~ n is defined in the context of Eq.(^). This choice guaran- 
tees that the S- matrix calculated field-theoretically to first order in w(K) and calculated 
quantum-mechanically to first order in v(K + , K±) are identical. The S-matrix carries a 5- 
function for light-front energy K~ between initial and final states. The definition of Eq. ( |63| ) 
removes that 5-function from v(K + , Kj_) and allows for general off ^"-shell matrix elements. 
Thus, Eq.(|64"D implies a very particular off-shell extension. This procedure of identification 
- it is no derivation - is standard for the instant-form of quantum mechanics, e.g., when 
the one-boson exchange potential between nucleons is introduced. This paper extends that 
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procedure to light-front quantum mechanics. Furthermore, the potential is usually defined 
in the two-particle cm. system, i.e., for K± = 0, and is considered unchanged in moving 
systems, i.e., independent of K± and K + . 

The identification fl63|) motivates a quantum-mechanical potential. It does not attempt 
to derive it. The goal of the identification is to simulate exact solutions of the BSE in 
best accord with a chosen physics criterion, quantum mechanics description has different 
objectives than matching a field-theory result. It rather attempts to describe many-particle 
systems with the same rules once it has done so satisfactorily for the two-particle system 
with the same rules. Thus, when the quantum-mechanical potential cannot be derived 
completely, as is the case in hadronic physics, the potential is tuned to known experimen- 
tal properties of the two-particle system and then considered a vehicle which carries that 
two-particle information to many-particle systems. Despite the particular many-particle 
aspect of quantum mechanics, a study of its predictive quality even for the two-particle 
system is interesting. Figs. 2-4 perform such study for the two-boson system of Sect. IV. 
The bound state constitutes an especially stringent test. For the instantaneous choice, the 
approximation, K~ = K~ n > m\ + m 2 in the interaction in the cm. system is quite 
severe, because in this case field theory requires K~ = Kg < mi + m 2 . The relation 
between the coupling constant gs and the bound-state mass Mb and the dependence of 
the bound-state wave-function f{q) on the momentum q = are compared in the field- 
theoretic and quantum-mechanical descriptions. Results are studied for the approximations 
w(K) ~ w( 2 '(K) and w(K) ~ w^ 2 \K) +w^(K) up to second order and fourth-order in the 
coupling constant gs- The quantum-mechanical binding energy and wave-function preserve 
most field-theoretic characteristics, expectedly better in the case of small binding rather tha 
in the case of strong binding. The quantum-mechanical choice of the potential is usually 
based on the one-boson exchange, i.e., on the approximation w(K) ~ w^ 2 \K). We are 
happy to find that this identification accounts better for the field-theoretic results than the 
choice based on w(K) ~ w^(K) + w^(K). 

Instead of solving Eq . (|60|) , its formal identity with the energy eigenvalue problem for a 
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nonrelativistic hamiltonian is often exploited |13j and \<Pb) is applied directly in the frame- 
work of light-front quantum-mechanics. 

The response of the quantum-mechanical system to an electromagnetic probe is given by 
a four- vector current j£(K' + — K + , K' ± — K±) which, as the quantum- mechanical potential 
is a three-dimensional operator and it depends on the three-dimensional momentum transfer 
(Q + , Q±) = (K' + — K + , K' ± — K±). As in the case of the potential, contact can be attempted 
with the corresponding field theory. A possible identification is 

(k'+k' 1± f v {K' + - K+, K' ± - K ± ) k+h±) := 



with 



i Wi^W) {h ^ mK ~ K)l kt ^ kHK K+ -kt) (65) 

K v = (K on ,K + ,K ± ) 

K' v = (K' onl K' + ,K' ± ) . (66) 

The field-theoretic j^(K' — K) is the one of Eq.(pETD in Sect. III. It contains the field-theoretic 
interaction in the form of w(K). The quantum-mechanical current j£(K' + —K + , K' ± —K±), is 
derived in the special case of elastic scattering between bound states. Thus, the identification 
of Eq . (165|) is not consistent with the choice of Eq.(^), which guaranteed the agreement of 
the field-theoretic and the quantum-mechanical S-matrix in first order in the interaction. 
Nevertheless, the quantum-mechanical current j£(K' + — K + , K' ± — K±) can be meaningfully 
studied and separated into interaction-free single particle and interaction-dependent two- 
particle pieces. Thus, the definition of Eq.(^) implicitly contains a possible quantum- 
mechanical definition of an interaction-dependent two-particle current. At this stage, the 



standard definition based on the identification of the S-matrix could also be given [ H| . It also 
exploits the formal identity of the eigenvalue problem with a nonrelativistic hamiltonian with 
equations like ((60), but it then identifies the nonrelativistic bound state with the solution 



b) of Eq.(B^p, the auxiliary field-theoretic bound state for the BSE. Thus, the calculation 



of the electromagnetic deuteron form factors in Ref. |T4| is performed in the field-theoretic 

19 



spirit of Eg. ((59|) . The two-particle current operators of pion range in Ref. JT4| should not be 
confused with the quantum-mechanical interaction-dependent two-particle currents of this 
section. 



VII. CONCLUSION 

The paper suggests a calculational procedure for solving the BSE with comparative ease 
and in principle, with any desired accuracy. The procedure is based on an auxiliary three- 
dimensional integral equation, in the framework of light-front dynamics, whose solution then 
yields the result of the BSE by quadrature. The intermediate auxiliary quantities do not 
display covariance; covariance is restored in the final step to the full result of BSE. 

The calculational procedure is exact, but it also offers an efficient approximative scheme: 
Only particles propagate. Antiparticles do not. Antiparticle propagation is relegated to 
the effective interaction. The convergence with repect to the number of exchanged particles 
mediating the interaction appears to be rapid. Though only an indication of that fact comes 
from the simple test case of a BSE bound state in ladder approximation, it gets supported 



by the similar result of Ref. |12| for the corresponding scattering amplitude. Calculational 
improvements are possible in a systematic manner. Thus, as a further and physically more 
interesting consequence, the solution of the BSE for bound state and scattering up to fourth 
order in the coupling constant, i.e., in ladder and crossed ladder approximation and with 
the inclusion of self-energy corrections is obtained based on a simplifying three-dimensional 
calculational procedure. The procedure capitalizes on beneficial properties of light-front 
dynamics. It should be an interesting alternative to the Gross approach || which is also 
three-dimensional and which has been suggested to include the cross-ladder exchanges ap- 
proximately. 

The calculational procedure is general, though it is given in this paper for an interacting 
two-boson system only. The ideas needed for an extension to fermions are developed but 
important technical details have not yet been worked out and unforeseen difficulties may still 
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arise. The problem of rotational invariance in light-front dynamics will become especially 
acute for fermions when spin and orbital angular momentum are to be coupled. The auxiliary 
three-dimensional quantities will then be hampered by their lack of rotational invariance. 
We strongly believe however, that the final step to the covariant result of BSE will overcome 
that difficulty 

The auxiliary thee-dimensional quantities, i.e., the operators and equations, that medi- 
ate the solution of the BSE, are close in spirit to relativistic quantum mechanics. The paper 
also discusses this relation. First, only particles and not antiparticles, propagate in the 
three-dimensional equations and in quantum mechanics. Second, the quantum-mechanical 
interaction is an instantaneous potential, the corresponding interaction w(K) in the three- 
dimensional equation is not. However, this paper finds that the instantaneous choice for 
the potential does not distort the physics of the underlying field theory. Thus, the relation 
between quantum mechanics and field theory can be made close. However, compared to 
field theory, quantum mechanics has the virtue of an instant extension to many-particle 
systems: Barring small corrections due to many-particle forces and the quantum mechanical 
interaction is additive in the instantaneous pairwise potentials. In fact, the conceptual strat- 
egy of quantum mechanics often is to tune away shortcomings of the chosen instantaneous 
potential by adjusting undetermined phenomenological parameters to vital known experi- 
mental properties of the considered two-particle system. In this way the potential carries 
the accepted knowledge on the two-particle system over to many-particle systems. 

The paper left open the relationship of the theoretical apparatus developed to realistic 
physics problems. We have in mind applications to hadronic and subhadronic systems. 
The concept of light-front wave-functions was applied in the context of nuclear physics to 



describe the deuteron [|T^] and the discussion of its properties in the light-front continues to 
the present |T5|| . The BSE is supposed to yield bound states and the scattering amplitude for 
those two-particle systems. In contrast, the response of such a two-particle system towards 
an eletroweak probe is considered in perturbation theory. The required matrix element 
is determined by the field-theoretic current between states of the BSE. The paper offers 
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two equivalent routes for calculation: Either the covariant states of BSE are constructed 
and then used in their four-dimensional forms or the field-theoretic current is reduced to 
a three-dimensional one, consistent with the auxiliary three-dimensional ones. Again, the 
latter calculational scheme is close to the quantum-mechanical one in spirit. The definition 
of two-particle exchange currents for the use in quantum mechanics is sketched. 
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APPENDIX A: EVALUATION OF AUXILIARY QUANTITIES 



The operators GQ{K)\gQ{K)~~ l and g^lGo^K) connect three-dimensional and four- 
dimensional basis states. The two operators are related by conjugation; we therefore discuss 
only one, i.e., G (K)\gQ(Ky l . 

The momentum space matrix elements of Go(K)\go(K) for K + > 0, are 

5 (k[ + - hi) S (k' 1± - k 1± ] 



y- _ "-ix" 



K- - K m - k? m + io) 0{K+ - kt )0{kj 

js- _ U- _ (K ± -k' 1± ) 2 +mj-io 
1 K+-k'+ 



X 



(Al) 



When the available light-front "energy" K is not on-shell, i.e., K ^ k lon + k 2on , the 
evaluation of the matrix element in Eq. (|A1 ) is standard. The two singular propagators 



k[ — fcix ^ 1 w j and (k — k[ — — ^+\y+' 2 J can be rewritten as a 5-function 



and principal-part singularity; integration on k'{~ can be carried out with usual techniques. 

A problem arises, when the available light-front "energy" K~ is on-shell, i.e., K~ = 
Kon = ^ion + ^2on- Without losing generality, we will have suppose that K + > and k^ > 0. 
Then, K~ — k{ on — A;^" on + io = +io and the limiting process of going to the real axis must be 
performed with care. However, in this situation the matrix element will always be integrated 
with respect to k'{~, over a function f(k'{~) still to be determined and, unfortunately with 
unknown analyticity properties, i.e., 

dk[-f(k[')(k[-k[ + k[ ± jGoWM^r 1 ! k t^) = r 5 ( k ' i+ ~ fc 5 - ki± ) 



2tt 

x / dk'r 7 V7 ^ (K- - k^ on - k 2on + io) . (A2) 

J 1 k'--k.7 4- in\ (K- - k'r - k- + V l0n 2 ° n ' K ' 



k[ - k lon + io) [K - k[ - k 2on + io) 

Without any loss of generality, we can think of f(k'{~) as being split into a part f u hp{k'i) 
having singularities only in the upper half fci~-plane and a part/i/ ip (/c / 1 ~) having singularities 
only in the lower half fc'^-plane, i.e., 

f(K-) = f uhp (k'r) + fin P (K~) • (A3) 
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In the case that there is a part with poles simultaneously in both half planes, they can be 
fully separated, i.e., 

1 1 



9( k 'i )tt 



k[ — a,\ — ia 2 k[ — fa + if3 2 
1 f 1 



(a-/3) + i(a 2 + P2) 



1 



(A4) 



k[ — ai — ia 2 k[ — Pi + ip 2 

with g{k'i) being singularity free. The integration in Eq . can now be carried out using 
Cauchy's theorem: 

J dk'ff(k'r)(k'rk'+k'i ± \G (K)\g (K)- l \ kfki ± ) = 5 (k'+ - kf) S (k' 1± - h 
— k lon — k 2on + 



LL X 



fuhp(k lon ) 



to) X 

1 



K k lon k 2on + io K k 2on k lon + io 



5 (k'+ - k+) 5 (k[ ± - k x± ) [f uhp (Kj + fih P (K 



h 2on/ 



(A5) 



We note that propagators cancel and no singularity remains. However, the result ( |A5| ) 
is for practical purposes useless, since the split into two parts with disjoint singularities 
is not known in a numerical calculation. If, however, the light-front "energy" is on-shell, 
K~ = K~ n , then the two terms can be recombined to the original function, i.e., 

J dk'.-fik'mk'rk'+k^ jGoWM^r 1 ! k ^ = 6 ( k ? - fc 6 - fQzJ ( A6 ) 



for K =K on . 
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APPENDIX B: INTERACTION IN FIRST ORDER 



The interaction w(k), defined by Eqs.(^) and (O) to lowest order of the driving term 
V(K), is given by 

w^\K) = g (K)' 1 \G (K)V(K)G (K)\g (K)~ 1 , (Bl) 

where the matrix element of the operator \G (K)V (K)G (K)\ is 

(k[ + k' 1± \ \G (K)V(K)G (K)\ \ kfk 1± ) = [ dk'fdk^ 

(2n) J 



1 1 1 



x 



k[+(K+ - k>+) ^- _ 3W-*>) [k- - k[~ 



(K±—k') 2 +m,2—io 



K+-k'+ 



1 



1 



^ + - k t) ( K - %^t») ( K - _ K _ 



(^ ± -fc 1± ) 2 +m|-tc 
fc+ / V" '"J- K+-k+ 



(B2) 



The double integration in k~ in Eq. flB"2]) is performed analytically using Cauchy's theorem 
and the condition K + > 0. The integration is nonzero for K + > k'^ > and K + > kf > 0. 
Two possibilities also appear for a forward propagation. For kt > k'i, a a is emitted by 
particle 1 and otherwise absorbed: 

i9(K+ - k'+)6(k'+) 



(k'+k' l± \ \G (K)V(K)G (K)\ \ kfk 1± ) = {ig £ 



k^(K + ki~) [K k lon k 2on + 



10 



x e(kj - k[ + ) i + e(k'+ - kf) 



(k^ kxj (K k lon k 2on k' aon + io^J {k^ kf^j [K k lon k 2on k aon + 



to 



x - «(^ -*»«(*»_ , (B3 ) 

k l (K + k 1 ) iK k lon k 2on + io) 



where the light-front "energies" of the intermediate states of the individual particles are 
given by 

kf ± + ml 



K lc 
k\ c 



k>+ 

kf ± + ml 
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k' 



2cm 



v 2on 



y- 
kZ 



(K± - k' 1± ) 2 + ml 
K+ - k[ + 

(K± ~ fci±) 2 + mj 
K+ - kf 

(k' 1± -k 1± f + f , 2 

kf - k'+ 
(k' 1± -k 1± ) 2 + fi 2 



(B4) 



k' + — k + 

The global three-particle propagator for 1, 2 and o appears in Eq. (|B3|) , in two cases: when 
a is either emitted or absorbed by particle 1. 

The matrix element (k'f r k' 1± w^ 2 \K) kfki±) is obtained from Eq.( P3|) by multiplying 
both sides by the matrix element of the operator g (K)~ l from Eq. fl25|) . 

(k[ + k' 1± y 2 \K)\kik 1± ) = 

6{kf - k'i ) i 



= (Ws) 

+ (igs) 2 
= {igs) 



(kf — k'fj (K — k' lon — k 2on — k' aon + to 

6{k' 1 + - kf) i 

(kf~ — kf^j (K~ — k lon — k 2on — k~ on + w 
2 6(kf - k'x) % 



(kf - k[ + 



K- 



y2 



2 (K ± -k 1± )2+m 2 2 



k'+ 



K+-k+ 



k+-k'+ 



+ 10 



+ {igs) 



2 ~~ Ef 



k' + - k + 



K- 



(K ± -k\ | ) 2 +m 2 



K+-k' 



fc'+-fc+ r 10 



(B5) 
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APPENDIX C: INTERACTION IN SECOND ORDER 



The interaction w(k), defined by Eqs. flHTD and (ITT]) to second order in the driving term 
V(K), is given by 

w(K) ~w {2 \K)+w {i) (K) (CI) 

where w@)(K) is given by Eq. QE5D and 

w^\K) = g G {K)- l \G Q {K)V{K)G Q {K)V{K)G G {K)\g Q {Kr l 

- goW-^GoWVWGaWVWGoWlgaiK)- 1 . (C2) 

The second term in Eq.(C2) corresponds to the iteration of the interaction w^ 2 \K) 



9o (K) ~ x | G (K)V(K) G (K)V(K) G {K) \g Q {K)~ l = 

w< 2) g {K)w® . (C3) 



The matrix element of the operator \G (K)V(K)G (K)V(K)G (K)\ is 
(k[+k[ ± \ \G (K)V(K)G (K)V(K)G (K)\ \ kfk 1± ) = J dk'fdp^dk^dp+d 2 p 1± 



1 1 

x 



k't (k+ - K + ) _ ga (k- _ K 



(K ± -k' 11 ) 2 +rn\-io 



K+-k'+ 



1 



1 1 



i i 



1 1 



(C4) 



kf (K+ - kf) ^- _ gx^iz!£ ^ _ fc- _ (£iz|ili^lz^^ 

The on-energy-shell values of the light-front minus momentum in Eq. (|U4|j are given in 
Eq.([B4|), and 
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Pic 



— >2 . 2 

Pi± + ™i 
+ > 
Pi 

(K± - Pi±) 2 + m\ 



P^on = - ~ K+ Z' pt * ■ ( C5 ) 
The matrix element (k?k' 1L \ \G (K)V(K)G (K)V(K)G (K)\ \ kfk 1± ) is found by ana- 
lytical integration in the light-front "energies" in Eq. (|C4|) . To separate the intermediate four 
particle propagation, which occurs for k't~ , pt and kf satisfying < kt < pt < k' + < K + , 
the following factorization is necessary 

1 1 



K- 



Pi 



k+-p! 



Pi ~h 



Cl-L— PI ) +H*—tO 



K- - K - 



K±-Pi±) +m|-io 



K+-p1 



X 



1 



K- 



Pi 



.Rx—Pix) +rrii—io 



+ 



1 



Pi ~h 



(C6) 



K+-p+ ri -l pt-kt 

After the Cauchy integration in the light-front "energies" the result for 
(k'-^k^ | \G (K)V(K)G (K)V(K)G (K)\ | kfk 1± ) in the region of < kf < pt < K + < 
K+, which is denoted by (k[ + k[ ± \\G (K)V(K)G (K)V(K)G (K)\ (a) | kfk 1± ), is given by 

(k'+k' 1± \G (K)V(K)G (K)V(K)G {K)\ [a) k+k 1± ) = 



2(2tt) 



dpfd 2 pi± 



6(k' 1 + )6(K+ - k'i) 



k[ + {K+ - fcj+ 



K- 



1*1 i 



Ki-kl 



K+-k 



+ io 



x [F'(K) + F"{K)\ 



6{kt)6{K + - kt 



k?(K+ - k+) 



K- 



K+-k+ 



(C7) 



+ io 



with 



F\K) 



6(k[ + -pj 



(K + - Pt 
6(pf)6(K + - pj 



K _ g^+mf (K±-k' 1± ) +m 2 (k' 1± -P!±) +P- 2 , ■ 



X 



pi(K+- P f 

6(pf - kf) 



jr_ _ pf_!_+mf _ (K ± -p 1± ) +m 2 



pT 



K+-pJ 



(pt ~ kt 



K- 



k 2 ± +m 2 (K±-Pi±) +m 2 +H 2 . 

kf K^Z-f + 10 
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F»(K) e{ ^-rt 



K+-k'+ K -Pi 

i 



jr_ kf ± +rrq (K^-k'^f+ml (k' 1± ~p 1± ) 2 +fi 2 (pix-fcix) 2 +fi 2 . 
K kf kT^t jr^T + 10 

x ^ ~ kl) i m 

(Pi - kf) k'f ± +ml (K ± -p 1± ) 2 +ml (p 1± -k 1± f+^ . 

K k+ KTZ^f -rz^ + ™ 

The part of the propagator given by Eq.( |<J7[ ) contains the virtual light-front propagation 
of intermediate states with up to four particles. The function F' contains only intermediate 
states up to three particles and is two-body reducible. It will eventually be canceled by 
the corresponding piece in the second term in (|C2|) . The function F" has one intermediate 
state in which four-particle propagator which can be recognized as the middle piece of 
Eq . ( |C9|) . The other possibility which includes up to four particles in the intermediate state 
propagation is given by < ky < p\ < kf < K + . To obtain this part, we perform the 
transformation k[ <-> k\ in Eq . (|C7|) . 

The contribution of the region < Pi < k~i < k'i < K + to 
the matrix element (k[ + k' 1± \ \Gq(K)V(K)G q (K)V(K)G q (K)\ \ kfk 1± ) is denoted by 
(k^k^ \Go(K)V(K)Go(K)V(K)G (K)\( b) kfhx)- It contains only up to three-particle 
intermediate states and is two-body reducible. Consequently it will be canceled by the 
corresponding piece of the second term in (|Cj|). It is given by 



(k[+k[ ± \Go(K)V(K)G (K)V(K)G (K)\ {b) k+k 1± ) = 

{igs) A f, + , 2 6(k' 1 + )6(K+ - k[ + ) i 

2(2*)* J Pl(1Pl± k' 1 + (K+-k[ + ) k? x+ mi (K ± - K± ) 2 + ml 

K fc^ K+^k? + W 

6{k'+ - fc+) i 

X (K + - k l)17- p1x+W (K ± -k> 1± ) 2 +m* (fc^-p 1± ) 2 +M 2 ■ 

6(pt)9(K + - pt) i 

x 



X 



pt{K+-pf) Pl ± +mf (K ± -p 1± ) 2 +ml 

O(kt-Pl) • 

{kt-Pl) K - (i?x-fcix) 2 +ml (fc!X-pix) 2 +M 2 , „■„ 

K K*=k? W^Pl +W 
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6(kf)6(K + - kt 



kf{K+ - kt 



K- - 



kj ± +mj (i?x-fcix) +r. 



(CIO) 



k+ K+-k+ 

For the momentum region satisfying < k'±~ < ki < pi < K + , the contribution 
to the matrix element (k?k' lL \ \G (K)V(K)G (K)V(K)G (K)\ \kfk 1± ) can be obtained 
from Eq. flClOl) by performing the following transformation on the kinematical momentum: 
k[ <-> K — k'i, k\ <r-> K — k\ and m\ <-> m-i. From Eqs. ( |C9|) and (|C10| ), the following result 
is obtained 



(k'+k' 



i± 



(k'+k' 



i± 



\G (K)V(K)G (K)V(K)G (K)\ \ kthj_) = 
\G (K)V(K)G (K)V(K)G (K)\ (a) \ A;+£ 1± > + [k[ <-> h 
+ ((k[ + k' 1L G (K)V(K)G (K)V(K)G (K)\ (b) \ kfk 1± ) 
+ [k[ K — k[, ki «-> K — ki, ni\ <-> m^) . 



(Gil) 



The subtraction of the iterated first order driving term in Eq. ( |C2| ) cancels the correspond- 
ing terms in Eq.( CJll ) such that the matrix element (k'+k^ w^ A \K) kfki±) is two-body 
irreducible with a global four-body propagation. It is obtained from Eqs.(|C7D, ( CT|) and 
(ED as 



(k'+k' l± w&{K) ktku) 



2(2tt) 



e(k[ + - P t 



(k[ + -pt) 



K- 



A"+-fc 



+ io 



x 



x 



_ fc ix+ m l 



(i?X-%j_) +r 



(fe' ix -pix) +M 2 (pix-feix) 



+ 10 



(pi - kt) 



K- - 



fc 2 ± +m 2 (K±-Pi±) +ml (pix-fcix) ■ 
~kf K+^p? pf=fcf + 10 



+ [k[ «-> fci] 



(C12) 
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APPENDIX D: INTEGRAL EQUATION FOR THE BOUND-STATE 



In the approximation considered, the vertex function satisfies an integral equation with 
the kernel containing two parts, one corresponding to Eq. (|B5|) and the other to Eq. (|C12j ). 
The plus momentum are rescaled by K + , such that the momentum fractions x = 

V — ~k+i ~ — k 

homogeneous integral equation for the light-front vertex function is evaluated in the center 
of mass system, 

1 r d 2 k 1± dx & 2 \y,P 1± ;x,ki ± )+K;( 4 \y,k' lx ;x,ki^ 



jkp, and z = jh:, are used. The notation (k'^k[j_ 7b) = 7s(y, k' l± ) is introduced. The 



7s(2/> k i 



id 



(2?r) 3 J 2x(l-x) 
where the free two-body mass is M { 



x(l—x) 



Ml-Ml 



and < x < 1. 



-1b{x, ku.) , (Dl) 



The part of the kernel which has only the propagation of virtual three particles states 
foward in the light-front time is obtained from Eq. (|B5|) , 



l-LJ 



gl- 



9(x - y) 



b' 2 

r2 K _i± 



1^ _ (k' 1± -k 1± ) 2 +H 2 

x-y 



+ 



X 



k 



(D2) 



(x - y) [Ml 

Eq.( pi| ) with the effective interaction given by ( |L)2| ) corresponds to the Weinberg equation 
derived from the BSE in the infinitum momentum frame ||16|| . It has also been solved in Ref. 
1 171 an d in Ref. [|0| including self-energy correction. The equivalent equation for fermions 



has been discussed in Ref. 19 



The contribution to the kernel from the virtual four-body propagation is obtained from 
Eq.flCm 

fc (4) (y,Ki,x,ki±) = 



d 2 p±dz 



(2vr) 3 J 2z(l -z)(z- x) (y - z) 

1 



9(z-y)9(x-z) 



Ml 



ft i± 



pf ± +rn 2 



(k' 1± -p 1± ) 2 +n 2 
z-y 



M 2 B -^ 



k 2 ± +m 2 



z~y 



(Pi±-fclJ_) 2 +/^ 
x—z 



M2 _g i± m ! 



k' 2 



1-x 



(pi±-k 1± ) 2 +p 2 



+ 



X 



y,k 



k[ 



(D3) 
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Eqs. (|Dl| )-( |D3|) are easily recognized to be covariant under kinematical light-front boosts. 
However, the covariance of the four-dimensional wave-function fl¥I|) is certainly lost by a 
finite expansion of W(K) in Eq.(|ll|) and the use of the corresponding w(K) while covariance 
continues to hold for the solution of Eq. (|TT|) . 
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FIGURE CAPTIONS 



Fig.l. Light-front time ordered diagrams for w^ 2 \K) (a) and w^(K) (b), representing 
the light-front time ordered view of one and two a exchanges, respectively. 

Fig. 2. Results for gs as a function of the two-body bound state mass Mb for fi = 0.5m. 
Numerical solution of the covariant four- dimensional BSE ([EJ) (solid curve), the light-front 
Eq.([3~8l) with interaction including up to three-particles in the intermediate states, i.e., with 
w(Kb) — w^(Kb) (dashed curve) and including up to four-particles in the intermediate 
states, i.e., with w(K B ) — w^(K B ) + w <yA \K B ) (dotted curve). Solution of the quantum 
mechanics squared mass eigenvalue equation (|60"D, with w(K v ) ~ w^ 2 \K v ) (long-dashed 
curve), and with w(K v ) ~ w^ 2 '(K v ) +w <yA \K v ) (short-dashed curve) defining the two-particle 
potential in Eq. (|63|) . 

Fig. 3. Results for the transverse momentum distribution f(q) as a function of the 
transverse component, q, of the individual four-momentum, for Mb = and fi = 0.5m: 
(a) numerical solution of the four- dimensional BSE with g s = 20.14; (b) relative error 
of the various approximations with respect to the four-dimensional BSE results, defined 
by Df(q) = 1 — f$ p {q) / ' fexactW) with n = 2 and 4. Results for the light-front Eq.(^8|) 
with an interaction including up to three-particles in the intermediate states, i.e., with 
w(Kb) — w^ 2 \Kb) where g s = 20.8 (dashed curve) and with an interaction including up 
to four-particles in the intermediate states, i.e., with w(Kb) — w^(Kb) + w^(Kb) where 
g s = 20.2 (dotted curve). Solutions of the quantum mechanics squared mass eigenvalue 
equation (|60|) , with the two-particle potential in Eq.(|63"D defined by w(K v ) ~ w^ 2 \K v ) 
where g s = 15.7 (long-dashed curve), and with w(K v ) ~ w^ 2 \K v ) +w <yA \K v ) where g s = 14.9 
(short-dashed curve). 

Fig. 4. Results for the transverse momentum distribution f(q) as a function the of 
transverse component, q, of the individual four-momentum, for Mb = 1.98m and /i = 0.5m: 
(a) numerical solution of the four-dimensional BSE with g s = 9.03 ; (b) relative error of the 
various approximations in respect to the four-dimensional BSE results, defined by Df(q) = 
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1 ~~ fappio) I ' fexactio) with n = 2 and 4. Results for the light-front Eq.(^) with interaction 
including up to three-particles in the intermediate states, i.e., with w(Kb) — w^ 2 \Kb) 
where g s = 9.10 (dashed curve) and with an interaction including up to four-particles in 
the intermediate states, i.e., with w(Kb) — w^ 2 \Kb) + w^(Kb) where g s = 9.03 (dotted 
curve). Solutions of the quantum mechanics squared mass eigenvalue equation (0), with the 
two-particle potential in Eq.fl63|) defined by w(K v ) ~ w^ 2 \K v ) where g s = 8.33 (long-dashed 
curve), and with w(K v ) ~ w^(K v ) + w^\K v ) where g s = 8.23 (short-dashed curve). 
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FIGURES 



Fig. la 



Fig. lb 
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FIG. 2 
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FIG. 3a 
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FIG. 3b 
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FIG. 4a 
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FIG. 4b 
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